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INTRODUCTION 


























| ۱ ar The use of mathematical models in energy system planning has become 
y Planning Studies utilize the 

i almost commonplace in the past few years, and both optimization and simulation 
essment or policy analysis, 0 3 
and regional 


analyses of global,” national” 


Reliability Council Region level. energy-economic systems have 


found widespread application in policy studies. Although such energy system 


utility to State Agencies, or to 


resource constraint questions models undoubtedly afford valuable insights to improved public policy making, 


ocedure to translate rezional especially where linked with Input-Output and/or econometric models of the 


h emphasis on the electric and M the difficulty of integration of environmental impact considerations 











considerations as an inherent Thus far the typical approach is one based on a computation of aggregate 
residuals; energy-economic models are used to determine the optimal technology 
mix; emission coefficients are applied to each system component for pollutants 


of interest (S0,, particulates, etc.); and the environmental impact expressed 







as total emissions. Yet it seems clear that such numbers are of questionable 






value for environmental policy analysis, let alone for any credible examination 








of environmental-energy-economic trade-offs. Attempts to refine such computations 






by replacing the usual objective function specification of cost minimization 


















by minimization of emissions (weighted, say, according to health impact and 





constrained by some investment limit) must also be regarded as somewhat dubious 





7 
surrogates for proper environmental analysis. 




















If one examines health impacts as a more specific example, the difficultv 
















is readily seen to lie in the fact that health impacts are largely determined 
by the exposure of populations to ambieat levels of pollutant concentration, 
not to aggregated emission levels: and ambient atmospheric (or water) concen- 
trations can be determined only by application of appropriate environmental 
transport models to specific point source emission data. The essential feature 


of such models, of course, is their locational specificity: a computation of 








atmospheric sulfate levels, for example, requires knowledge of stack locations 





to a considerably more exact level than the level of disaggregation of currently 






available energy system-economic models. Thus the prediction that so many Gw 








of coal fired electric generation capacity with a given emission control tech- 














۱ ; 8 | i 
nology will be located in, say, a particular Census or Reliability Council 


(typical outputs of multi-regional energy system nodels) still yields 















relatively little of policy relevant value. To be sure, the regional scale 
imparts more information than a national level aggregation, but one is still 


some way from che desirable information content. 






The importance of moving beyond mere analyses of residuals is perhaps 
nowhere as striking as with the key coal related health impact issues of 
sulfates and ozone. Both pollutants are known to be of significant health 
impact, and known to be correlated with the presence of coal combustion, yet 
neither substance occurs in significant amounts as a stack emission. Rather, 
local ambient levels are determined by complex atmospheric conversion processes 
that can be captured only by the application of regional scale meteorological 
models. Thus, sulfur dioxide (ر50)‎ and nitrogen oxide (NO) stack 
emissions are not amenable to simple linear transformation to sulfate (SO, ) 
and ozone concentrations, which in turn underscores the very limited value of 


50, and NO, residuals as measures of environmental damage. 


It is clear, then, that there is a need for bridging the gap between 
energy system scenarios on the one hand, and environmental transport models 
on the other; and this bridge requires a siting of polluticn emitting facilities 
to a level of geographic detail consonant with the spatial resolution of 
environmental models. It appears, moreover, that the level of spatial resolution 
that best represents the compromise between the ability to disaggregate a 
regional scenario on the one hand, and the necessity for locational specificity 
of environmental models on the other, is the county. This paper develops a 
linear programming model that achieves this link; in our case the specific 
design criteria are the ability to interface with the Brookhaven integrated 
multi-regional energy and interindustry aodel which determines the regional 
configuration of the electric sector, and the Brookhaven Diffusion-Irajectory 


E 


Model for study of regional air quality and health impact issues. 


It should be stressed that the requirements for a siting analysis 
addressing the rather distant time horizon of 2000 and beyond, the typical 
planning horizon of scenario analyses, are clearly different to the require- 
ments for a siting analysis within the normal planning herizon of electric 


utilities whose objective is to actually icentiiy specific sit 





= 
case, the rationale for siting facilities has little to do with assisting 










BASIC MODEL FORMULATION 
utilities identify future sites; 








































is dictated by the requ 









that for meaningful environmental analvsis of regional scale problems to be ا‎ TEC i | | 
On the assumption that the generation mix in the electric sector is 
conducted, facilities must be located in space at about the county level 2 


specified by the output of a regional scenario model, the siting problem can de 





The overriding objective is to demonstrate the relationship of 





























formulated as the determination of that set of facility locations that nininizes 
siting policy to environmental impact, and, as such, a siting exercise of 
the sum of resource extraction, transportation, and electric transmission costs, 
the type implied here will always involve at least two siting scenarios. 
subject to meeting certain energy and resource mass balance conditions. 
Interest is thus focussed not on whether, say, mine-mouth siting is "better" 
simplicity of exposition we shall assume here that only the extraction and 
or "worse" than load center siting of coal electric plants, (a question that 


transportation costs of coal are relevant--allowing, for example, an examination 

















would require a multi-objective analysis on a case by case basis, clearly ۱ 
of the trade-offs between mine-mouth and load center siting of coal burning 


beyond the logistic capabilities most types of assessment exercise), but 


1 , power plants. 
to illustrate what typical siting patterns would be if one were to follow 









either a mine-mouth or a load center strategy, and what regional scale impacts Such a problem definition proves to be readily formulated in terns 





would follow. of a linear program. In structure such an LP is partly a transportation model 


(with the possibility of fuel shipments from each producing untv to each 


consuming county), part transshipment model (with shipments of electricity only 









between adjacent counties). Although the transportation formulation for fuels 
considerably enlarges the size of the. problem, it is necessary to account for 
individual county-to-county shipments because fuels fre 
may have different heat, sulfur and ash content (whereas in electric transmission 


mo such distinction need be made). 


2.1 Generation Mix Constraints 
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denote the capacity of type x located in county 


is defined 
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Pumped Hydro 


nen, in order zo meet a regional (col : we requi 
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ach coal technology is related to the capacity bv the relationship 
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(e) 8760 5, (2) 10° Ts‏ 5 لس 
)13( یں + یں w Pop Farer‏ ہہ 
i (hrs Btu || KW‏ 
T mmm‏ - 
GW [year Kuh | | GW‏ 
K‏ © ۶ 
h; ۳ = 0 tel,...T for all 1 (14)‏ 


eration aix equation for coal (Eq. 2) aust be revised as follows 


cilities 


the model to a real situation, cognizance must of course be 

capacity. Thus, for example, in a model run for a vear 2000 
art of the necessary capacity is already in place, or would 
say, 1985. In effect, chen, one must constrain some of the 
heir 1985 values (minus any planned or likely retirements 


2000). This is accomplished by adding a series of constraints 


k 8 یف‎ 7K " 
ې‎ “Y, “or any X, 20 (15) 


installed capacity of tvpe k in county i in the preceding 
.0 year 
د‎ capacity retirements of type k in county i in the interval 


ا 


1 the preceding scenario year and the year of the model run 


(00 model run, che XŠ, Y refer to 1985 and the 1935-2000 


-tivelv; and for the 2020 model دا‎ to 2000 and the interval 





For existing transmission facilities, the most obvious constraint 
the type 


AU. Y.. for an Ti: 
IQ œ +) 7 š 


- 
- 


where ود کپ‎ the transmission capacity between i and j, proves to de inappropriate, 


since if one examines the transmission trid map, there is no indication tor a given 
link (1,4) as to whether power flows from i to j or j to i. Inspection of the load 
flow diagram would of course give an indication of the flow direction, but since 

the load flow analysis is available only for the existing year, imposition of 
additional transmission facilities (as calculated by the model for the scenario year) 
makes it quite likely that significant load flow changes will have occurred. Thus, 

a flow from i to j in the existing year by no means assures a sinilar flow direction 


in the scenario year. 
Fortunately, however, close inspection of the model formulation reveals a 


simple method for resolving this difficulty. In fact for a given node pair (3,37: 


only one of the two possible ¿low directions can be non-zero; thus 1f 


Y » 09, ۲. > 01 06 SET > ©, ود‎ > O. Non-zero values of both Y,, and 
1j ji ji ij +J 


Yi, would imply a closed loop, which cannot occur. Thus in place of Eq. (16), we 


may write 


y 7 
Fay ^ w (17) 


since the slack ín this equation is the in capacity required (in che event 


` 


that Ti, is inadequate), this can de writt 


the scenario 


remains only to determine 
existing vear transmission grid. Inspection readily 


circuits í a county line. Then, armed with t: lari impedance ioa 
curve and the cista 3,, Jerween 
SI 


between tne 












The Objective Function 
















Having introduced the problem variables, we are now in a position to 


State the objective function to the basic transhipment problem. If 





Yi, = extraction cost of 2-th coal in county i 







extraction cost of 2-th coal in the p-th province 









Y = haul cost from province p for tha i-th coal, $/ton/mile 







the set of gateway counties for province p 



























۵ = capital cost of transmission, $/GW/mile 








Y m coal transportation cost, S/ton/mile 
= transmission loss factor, miles ^l 


وو و 
heat content of the 2-th coal, tons/i0 btu‏ = 8 
cost per kWh lost in transmission‏ 
















CRF 


applicable capital recovery factor 






then the objective function can be stated as 













(N) $ 


Minimize 5 = CRF ۵ Y ka 
c ij 1j 


۶ 5 
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Annualized Capital Costs of New Transmission Line Requirements 
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Cost of Coal Shipments Within Cost of Imported Coal 
The Region 


3. ENVIRONMENTAL CONSTRAINTS 


3.1 Air Qualitv Constraints 













Air quality requirements are among the most important determinants of 
fossil plant siting. Two types of constraints can be imposed on the model. 
The first takes into account New Source Performance Standards (NSPS), which 
are statutory limitations on the emissions from new facilities." Because 
many State Implementation Plans for meeting Clean Air Act requirements 
specify NSPS at the county level that may te lifferent (and more stringent) 
than the Federal Standard (as shown on Table 1 for the case of so, standards) 
spatial variations in NSPS have bearing on the location problem. The second 
category of constraints captures ambient air quality requirements, which re- 
quires that emissions be converted to ambient concentrations using appropri- 
ate atmospheric modeling tools. The present version of the siting model 
includes only consideration of 50. emissions, and this section develops the 
constraint equations for these pollutants; however, the addition of particu- 
lates, NO, or hydrocarbons would follow exactly analogous formulations. 

Incorporation of the New Source Performance Standard proves to be 
relatively straightforward. Defining 


Xe ° so, Removal Performance of Subtechnology t 


8, = Heat content of the ¿-th coal, in tons/10”” Btu 


g, ٭‎ Sulfur content of the 2-th coal in lbs/lb coal 






Theá SO. emissions in county i, in 15 SO,/year, are given by 


2 


r 1 2 . جم‎ ۰ 20 ۰ $ L سه‎ |» 
“ir "1 "1 2000 398 ( Yo (20) 


ري« 





But the SO, emissions will be limited by the new source performance standard 
{ 6 


for that county, say ¢, 1b so,/10 Btu. Hence 





















Connecticut 
Massachussetts 


New York 


New Jersey 


Pennsylvania 










District of Columbia 


Counties in the Northeast with New Source 


Performance Standards for SO 
Fueled Electric Utility Boilers below the, 
Current Federal Standard of 1.2 15 50 0 Btu 


County 


All Counties 


Norfolk 
Suffolk 
All Others 
Bronx 
Kings 
Nassau 
New York 
Queens 
Richmond 
Rockland 
Westchester 
Bergen 
Burlington 
Camden 
Essex 
Gloucester 
Hudson 
Mercer 
Middlesex 
Monmouth 
Morris 
Passaic 
Salem 
Somerset 
Union 
Allegheny 
Beaver 
Bucks 
Chester 
Delaware 
Lawrence 
Montgomery 


Philadelphia 


from Fossil 

















lb 50,/10° Btu 
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One such constraint will axis: Sar each isc 


feasible. 


The incorporation of ambient air qual: 
and requires prior knowledge of 2 
characteristics, and explicit, quantitativi 
matrix, say T, che element T of which de: 
concentration in county j attributable to 
Ceteribus Paribus, Ti, will be a function < 


stack, the greater the geographical disper: 


Recall the expression for total SO, emissi 
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The determination of the inter-county transfer matrices رر‎ is not a minor 


matter, since it requires application of a trajectory-diffusion model to 

each county in the problem area. Because the computation of sulfate concen- 
trations depends on complex atmospheric chemistry, and because such 
diffusion-trajectory models require the solution of three-dimensional partial 
differential equations over long time periods (to simulate a typical month's 
or season's weather), the computational effort is considerable. However, 

once these calculations have been completed, and the transformation matrices 
for a unit emission determined, the matrices can be stored on tape and accessed 


by the matrix generator for the LP as sata.” 






3.2 Coal Cleaning 





Coal cleaning, designed to remove by mechanical means some of the 
inorganically bound sulfur (principally Iron Pyrites, FeS), is beginning to be 
used on a quite widespread basis as an alternative or adjunct to fiue gas 
desulfurization (FGD). Indeed, there are many coals that can meet current 
State Implementation Plan requirements for sulfur emissions by coal washing 
alone, without a FGD system, and the issue of the degree to which coal 
cleaning qualifies as Best Available Control Technology under the 1977 Clean 


Air Act Amendments is now of some importance. 


Coal cleaning can be captured by the model by allowing a conversion of 
coal of type x to coal of type y. Thus, for a location i, we adjust the 


7 
coal mass balance equations by the quantity of coal cleaned, Ci as follows; 


C Pb T € e E € e (٢ C - C, = 0 (24) 
t ixt 4X 
ix jal iix jel ijx tel i 











Coal to Cleaning 
Plant 


and for ۷ : 
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Coal Output From 
Cleaning Plant 
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3.3 Water Resource Constraints 










where E? is the Btu yield of the coal cleaning process. To the 0 


function is added the term 











L ¿ C 


"MN 
ieC x,ye(X,Y) ix ce 





where C is the set of locations at which coal cleaning is permitted. 
(X,Y) is the set of coal transitions allowed. 


is the unit cost of coal cleaning, in 1069/1012 ۷ ۰ 














For application of the model in the Northeast, we allow coal cleaning of 
high and medium sulfur coals to produce coals of medium and low sulfur, 


respectively (i.e., we define the set (X,Y) as (1,2) and (2,3)). 



















where L is the number of river basins in the study region. 



















Consumptive use restrictions are of the more important siting criteria in 


inland waterways now that once-through is being replaced by evaporative 
cooling methods. These restrictions are best handled on a river basin basis. 
Suppose S, is the set of counties in the 2-th river basis, and let Uy be the 
be the consumptive use per Gw of nuclear capacity, in cfs/Gw, and Ur the 
consumptive use per Gw of fossil capacity, also in cfs/Gw. Then if the total 
maximum allowable consumptive use for power generation is u, cfs, we have 


constraints of the form: 





















Such consumptive use constraints may, of course, be too restrictive, 


since several methods can be used to get around such a constraint. One is by 
low flow augmentation, which requries construction of a resevoir. Another 
is by the use of dry cooling (assuming, of course, that once-through cooling 


in inland and estuarine areas continues to be precluded). 









Notice that this allows a fairly accurate simulation of reality. 


example, under current conditions and policy, there 





Introducing She new variables: Delaware River Basin, 












can be no consumptive withdrawals without offsetting storage, and thus we 








Low flow augmentation reservoir capacity, in 1000 acre-ft, 


۳ 
would simply set u. = 0. Also, if one wanted to examine the impact of noe 
5 


in the i-th river basin. 





allowing low-flow augmentation reservoirs in a particular river basin (or 


basins), one would simply delete the appropriate A, vàriébta(s)^". Alternatively, 

















Gw of dry cooled nuclear and fossil, respectively, in the 







L-th ° 
th river basin one can limit the amount of low-flow augmentation by addition of constraints 
















1000 Acre foot of storage required to develop an additional of the form 





cfs of low flow, in the 2-th river vedra!” 













Consumptive Use in Estuaries: Special conditions may apply to thermal 









Then (26) may be rewritten as 


plants in estuaries in respect to consumptive water use. In particular, as 






the location of a plant moves further toward the ocean, the amount of low flow 
augmentation to be provided upstream becomes less. For example, the Hope Creek 
nuclear plant, located in Salem County about midway between the start of the 
Delaware estuary (at Trenton) and the Ocean, has been required by the Delaware 


River Basins Commission to provide low flow augmentation in the amount of 18% 










of what it would have had to provide if it had been located further upstream 


24 
near Trenton. 











Eq. (27) can readily be modified to account for this special condition. 


If we define the set of estuarine counties in the 2-th river basin as E, > 








Since both the low flow augmentation and dry cooling involve incremental then (27) can be reformulated as 


costs over once-through cooling, the objective function requires the 


addition of two further terms, namely 








+ 
Dago) 















Cost of Reservoir Penalty for Dry 


A is the reduction in requisite flow augmentation applicable to the 








is the cost per unit of reservoir ca y 7 ۱ i i 
P capacity, in $/1000 AF (for i-th countv; hence the product consumptive use of freshwater 


which estimates are available, by river basin, from U.S. 


in county i. for Salem county would thus be .13. 






of Engineers Studies 4 ١ 
e define then we may rewrite (31) as 


is the cost of drv cooling at nuclear stations (see below). 












is the cost of dry cooling at fossil stations (see below). 

























for a plant of type k, 





design capacity if cooled by wet or dry towers, say D 
OTHER CONSTRAINTS 








3ut, further, the operation of the towers themselves 







cooled by towers of type o. 






requires diversion of some of that design capacity, say P 


Exclusion Constraints 


A+ 









net plant send-out, available to meet the electric demands of consumers, is only 























Exclusion constraints can be handled very simply in the model by deletion 






of those generation capacity variables that correspond to excluded locations. 

































Penalties for Evaporative and Dry Cooling Systems: Both evaporatively cooled 


and dry-cooled thermal facilities impose significant additional costs to elec- 


If the net plant capacity for a nominally sized 1200 Mw plant of type k, For example, for nuclear siting one might exclude all counties that lie 





cooled by towers 0 is No, k' “hen for l Gw of net output the nominal plant size is totally within a given seismic zone; if S is the set of excluded counties, then 













tric generation. To the extent that the model allows an examination of the 
we simply delete all 





trade-off between cooling mode and location, any cost penalty over once-through 









cooling should also be captured in the objective function. ^" Once-through 


penalty in KW, is 


ju 


per KW of capacity of type k is 


the capacity 
cooling, although currently precluded in Inland waters, is an allowable option / variables that are then no longer eligible for entry into the optimal solution. 






2 4 c Ä 14 4 , , . , 
for off-shore and possibly even some coastal facilities. Similar formulations can be used for population density exclusions, or land use 


Because consumptive use constraints are handled on a river basin basis, exclusions due to prohibitory legislation, Federal or State reserved lands 










(State Parks, National Forests, etc.), or for reasons of an ecological nature 


we shall introduce the following new variables 









If the capital cost then the annual penalty 





such as the presence of an endangered species. 















cost is given by 





"W Gw of Evaporatively cooled nuclear capacity in the 2-ch river basin 





4.2 Equity Constraints 











O 
"dq 









= GW of Evaporatively cooled fossil capacity in the 2-th river basin 


With increasing concern by State Regulatory Agencies over 6 justification 










= Gw of Dry cooled nuclear capacity in the i-th river basin 
n for additional capacity increments, and over the equity issues raised by large 










Also if the capital for cooling system construction of type 0 at a 


generating facilities located in remote areas serving distant, and often out-of- 


= Gw of Dry cooled fossil capacity ín the 2-th river basin 





facility of type k is given by 7 تفن‎ of design capacity, then 


state load centers, the ability to simulate mandated levels of energy autarchy is 







Gw of net output 







of some policy relevance. The most immediate constraint along these lines is 


penalty coefficients: 









to require that a certain fraction of the peak load in each state must be met 







Penalty, in 3/vr/Gw of evaporative 





This is readily formulated oy the expression 





by capacity within that state. 





















Penalty, in $/yr/Gw of evaporative 





















Penalty, in $/yr/Gw of dry nuclear 








Penalty, in $/yr/Gw of dry fossil 

























Hence, the total penalty cost is the summation of (33) and (34), namely 


fraction of each state's peak load to be satisfied by generation 


A number of factors must be considered in deriving the appropriate penalty 





capacity within the state 


First, a plant of, say, 1200 Mw cooled dy once-through would have lower 6 1200 7 
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Finally, we must modify constraints (2 souneses com 
slack variables 92: se two inequalities equal the evaporative 


P and 623) and (29) 






variables 






This constr: 
of control total: 
additions in eact 
area, although cc 
service areas. 
4.3 Other Refine 


The above e: 





an illustration « 
many refinements 
policy issues cat 
"counties" with : 
underwater cable 
possibilities; a 
of the assumptio: 
maximum capacity 
trade-offs, and 
more sophisticat 
sites; and incor 
istic variables. 
Finally, th 
lation and what 
most important i 
which sets it ap 
analysis of alte 
of Masse and Gib 
treatment of wat 
tions of Provenz 
focused on gener 
focused narrowly 


broad treatment 


5. Comoutation: 





The model 3 


int set could also be used for maintaining other types 
|; for example, it could de used to ensure that capacity that State, and the Pennsvlvania-Jersey-Marvland Interconnection (PJM), cov- 
utility service area be related to the peax load in that ering Pennsylvania, New Jersey, Maryland, Delaware and the District of Colum- 
unty boundaries may de difficult to align with utility f bia. Table 2 indicates some of the computational features for typical model 
runs; using the APEX -III LP package on CDC 7600 equipment results in a com- 
ments and Extensions putational burden that even at commercial rates is seen to be most reasonable. 
position elaborates only the most important constraints as These models have been used in a number of recent energy systems planning 
jf the methodological approach. As model development progressed, ۱ studies, including an analysis for the National Academy of Science Committee on 
and additional constraint sets were incorporated as specific Nuclear and Alternative Energy Systems (CONAES) as a basis for computations of 
ne under study. These included the addition of wash : air quality and health impacts, and the generation of siting scenarios for the 
appropriate special conditions for electric transmission (by ۱ National Coal Utilization Assessment, an ERDA sponsored effort to investigate 
| and exclusion sets to analyze various offshore siting ۱ : environmental impacts of increased levels of coal utilization. 
reformulation of electric transmission allowing relaxation ۱ ۱ As might well be anticipated, the major Problem in implementing such 
of constraint transmission load factors; constraints on : nodels lies in data assemblv, the development of the matrix generator that 
to allow examination of energy center and dispersed siting | i prepares the LP in a form suitable for the commercial LP code, and the in- 
to limit gas turbine capacity to VM HMM, 6 که‎ terpretation of the output-currently still relving on a somewhat tedious 
M population uy ses and their impact on nuclear ۱ manual transfer of the optimal basis to maps, but shortly to be automated 
poration of reliability issues by the addition of pea" | using 2 CDC developed matrix nanipilation code (MAGEN) and a TEKTRONIX 
| graphics terminal with harc copy unit, which will allow Maps (of the type 
ere are a number of key differences between this model formu- ۱ shown in Figures 1 and 2) to be directly generated by the LP. 
appear to ba quite similar Models an the literature. ihe ۱ As an exemplar of the type of siting scenario that emerges from use of 
s the fact that the generation mix is exogenously سیب‎ | T the model, consider the maps of Figures l and 2, which show siting scenarios 
art from the lineage of linear programming models "e ۳ for the New York Power Pool for the vear 2000. Perhaps the most remarkable 


rnative capacity expansion plans based on the pioneering i ; ١ ; è è ۱ ۱ w 
ture of these scenarios is the lack of fragmentation, since tae LP 


27 ۰ a . 28 ۲ ` ; 
rat” , and reviewed more recently in Anderson . | kasé ' ) tal ند‎ i 
EN lation captures neither scale economies or any =: charg a scatter 


sr resources it is perhaps more similar to recent 4 1 di 
29 X of small plants throughout the State might | : . However, even 
ano“ and Brill, et al” ; the former, however, is RES ۱ | | ند و‎ f l 

| the Figures do incorporate rounding up and zo + LOCMw, zane computati 


— < 


ation mix optimization, with little spatial detail; the latter is ae I ne ۱ l 9 I 1 
۱ experience thus far indicates that if a location is cost er-ective Ior à 
on water resource allocation problems, rather than on a more | ۱ : ۱ 1 | 
ticular capacity type, then that location will be filled to its maximum 
of energy facility siting. Lo 0 | l : ۱ i ۱ 
possible extent. A complete discussion ol these scenarios [rom an 7 
31 


and siting strategy point of view is to | found alsewnera . 
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Asvects 


s described above was implemented for chree power 3 ols 


0 
| United States: The New England Power Pool (NEPOOL), cover- 


England States, The New York Power Pool (NYP?) cover 
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1. See e.g. M. Kennedy "An Economic Model of the "orld Oil Market” Bell 
Journal of Economics and Management Science, 5,2 p. 540 (Autumn 74). 


2. See e.g. A. Manne "ETA: A Model for Energy Technology Assessment" Bell 
Journal of Economics, 7,2 p. 379 (Autumn 76) or E. Hudson and D. Jorgenson 
"U.S. Energy Policy and Economic Growth" Bell Journal of Economics and 


Management Science, 5,2 p. 461 (Autumn 74). 




















3. See e.g. Batelle Pacific Northwest Laboratories, "Regional Analysis of 15. 


the U.S. Electric Power Industry" BNWL-B-415, July 1975 or P.L. Joskow 
and M.L. Baughman "The Future of the U.S. Nuclear Energy Industry" 
Bell Journal of Economics, Vol. 7 No. l, (Spring 1976). Both these 
studies used Reliability Council Regions as the spatial unit. 


&. For discussion of the BNL-University of Illinois Model, see e.g. ۰ 
























Tessmer et al, "Coupled Energy System--Economic Models and Strategic 16. 
Planning" Comput and Ops. Res., Vol. 2, ۰ 213 (1975). 

5。 See especially P.F. Palmedo "The Incorporation of Environmental Consider- 
ations into Energy Policy Analysis" Presented at the AAAS meeting, Jan. 17. 
1975 (Published as BNL 19647, Brookhaven National Laboratory, Upton, 
New York 11967). ۳ 


















6. This is the basic philosophy of the Strategic Environmental Assessment 
System (SEAS) developed at EPA. 


7. See e.g: M. Beller "Energy Systems Studies Program" BNL 50539, Brook- 
haven National Laboratory, Upton, New York 11967, June 1976. 


8. Census regions are used in the Brookhaven Multiregional Energy Systems i 
Model. 





















9. Reliability Council Regions are used, for example, in the Baughman- 19. 


Joskow model (note 3, supra). 


10. R. Goettle, E.A. Cherniavsky and R.G. Tessmer "An Integrated Multi- 
Regional Energy and Interindustry Model of the United States", San 
Francisco TIMS/ORSA Joint National Meeting, May 1977. 












ll. R. Meyers and R. Cederwall "Fossil Pollutant Transport Model Development 
P. Palmedo, Editor, in "BNL Regional Energy Studies Program, Annual 
Report FY 1975", BNL 50478, November 1975. 

















12. The Electric Power Research Institute (EPRI), in a recent project on 
regional siting methdologies with the Applied Physics Laboratory of 
John Hopkins University, used not the term "siting" but Probable Dis- 

tribution of Effluent Sources from Energy Supply and Conversion”. 














Regional Science Association, 7th Annual Meeting, May, 1977 Halifax, N.S. 


Analysis for Environmental Management" International Regional Science | 
Review, Vol.l, No.l, Spring 1975, p.1-26 or, in an urban context, M. Pravin 


A Synthesis", Journal of Environmental Economics and Management, Vol.3, 


My colleague Michael Rowe has suggested a more refined formulation that 


stacks and other control devices. For three release levels, say corres- 
ponding to short, medium and tall stacks (h=1,2 and 3, respectively), 
Eq. (21) would be modified as 






20. 











Expressed differently, the trade-off between electric transmission and 
nuclear fuel transportation is overwhelmingly weighted to the former, 
implying that increased transmission distances necessarily increase 
total costs. 
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21. In the Delaware, for example, ? is about .420 10004F/cf3 - See e.g. 
U.S. Corps of Engineers, North Atlantic Regional Water Resources Study, 
Appendix R, Water Supply, p. R-12. 






















22. Ibid. 






Gateway countries are those through which coal is assumed to 5e imported; 
see Meier, note 15, infra. 



















23. In view of the demise of the Tocks Island Dam project in the Delaware 
River Basin, this would also be a not unreasonable subject for examination 


The assumption of a load factor constant for all transmission lines is, 
in the Northeast. 


of course, somewhat unrealistic. A more rigorous treatment of trans- 
mission is discussed in P. Meier, "A Linear Programming Model for 
County Level Energy Facility Siting", Brookhaven National Laboratory, 
Upton, N.Y., June 1977. The difficulty with a rigorous treatment is 
the significant increase in problem size. 




























24. See P. Meier and D. Morell "Issues in Clustered Nuclear Siting" BNL 50560, 
Sept. 1976, Chapter VII. 













25. E. is of course a subset of S.; the set cf non-estuarine counties in the 


P. Meier "Regional Science and Energy Policy: A Methodology for the i-th river basin is thus given by 5,- £,. 


Assessment of Increased Coal Utilization in the Northeast", Northeast 
































26. This method of approaching the cost penalties of closed cycle cooling 
systems in siting models was suggested by F. Moore in a series of unpub- 
lished memoranda to the BNL Regional Studies Program, available from the 
writer. 





They will vary from technology to technology; and are defined for a 
number of criteria pollutants such as so, particulates. 














27. P. Masse and R. Gibrat "Application of Linear Programming to Investments 
in the Electric Power Industry” Management Science, Vol.3,2, p.146-166 
(Jan. 1957). 


The notion of the inter-regional pollutant transfer matrix has been 
part of a number of theoretical models reported in the regional science 
literature; see e.g. W. Isard and R. Van Zele "Practical Regional Science 
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calls for a separate transfer matrix for each release height, a formula- 
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where Tan is the contribution of a unit emission at stack height h in 


county i to 50, concentration in j. This formulation, however, requires an 
additional cross classification for plant type, namely by stack height. 


The key assumption here, however, is one of linearity, and in particular 
in the case of sulfates, that the transformation of SO, to SO, involves 
first order chemical kinetic 2 " 
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